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The aim of this paper is to give a short overview for the set of publications consider- 
ing recently found effect of non-equivalence of multiplicity fluctuations in relativistic 
gases with globally conserved charge and energy. 

I. INTRODUCTION 

It was suggested to use the statistical approach to strong interactions more than 50 years 
ago pj, |2J, [3J. It appeared to be surprisingly successful in describing experimental results 
on hadron production properties in nuclear collisions at high energies (see e.g. Ref. 
{(J and references therein). This motivates a rapid development of statistical models and 
it raises new questions, previously not addressed in statistical physics. In particular, an 
applicability of the models formulated within various statistical ensembles. Recently, it was 
found that global conservation laws suppress multiplicity fluctuations and this suppression 
survive even in thermodynamic limit 7|, 8|. This unexpected result gave rise to the set of 
publications on this subject I7f-l22|. This paper gives a short overview starting from simple 

q n n n 

models [7|], [8j to the recently found experimental confirmation of this effect 19j, [201 ] . 



II. MULTIPLICITY FLUCTUATIONS 

Multiplicity fluctuations can be quantified by the scaled variance. For positively, and 
negatively, charged particles the scaled variance reads: 

, _ {Nj) - (N±f 

u = m ' () 

where angular brackets ( ) means averaging. The scaled variance is a useful measure, 
because for Poisson distribution it equals 1, independently of its mean value: 

^poisson 1 (^0 

Thus, the scaled variance says how much the studied system is different from Poisson distri- 
bution. Experimentally, the averaging in the Eq. ([T]) means the averaging on event-by-event 
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basis: a given observable is measured in each collision event and the fluctuations are eval- 
uated for the selected set of these events (see, e.g., review [(J). To calculate a statistical 
"background" for multiplicity fluctuations one has to choose a statistical ensemble for this 
calculation: grand canonical (GCE), canonical (CE), microcanonical (MCE) or grand mi- 
crocanonical (GMCE), see Fig. [U Usually authors do not make the difference between MCE 
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FIG. 1: Conservation laws in different statistical ensembles. 



and GMCE and call both as microcanonical ensemble. We introduce different names fol- 
lowing the suggestion of a referee, because we analyze and compare different ensembles in 
details. 

The choice of an ensemble depends on the experimental situation. If one exactly knows 
the energy, volume and charge of the system then such a system should be described in 
the MCE. Sometimes temperature of a system with exactly known electric charge can be 
measured much easier then its whole energy. Then such a system should be treated in CE, 
etc... In practice, calculations in CE and especially in GMCE and MCE are very difficult 
thus real calculations are always performed in GCE. One usually refers here to the textbook 
statement that all ensembles are equivalent in thermodynamic limit. 

This is the case for particle multiplicities. Different ensembles are equivalent if one 
choose a temperature and chemical potentials in a way that some exactly fixed variable in 
one ensemble equals to its adjoint average value in another ensemble, e.g. temperature T is 
defined from the condition E m c e = (E) c _ e _, and chemical potential /iq from the condition 
Qc.e. = (Q)g.c.e., etc., see Fig. CD However the equivalence of statistical ensembles does not 



apply to scaled variances. This was firstly found in [7| and will be illustrated below. 
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Canonical ensemble 



As a simplest example, let us consider a relativistic system in equilibrium which consists 
of one sort of positively, N + , and negatively charged particles 1 , N_, with total charge equal 
to Qc.e. = N + — iV_ . In the case of the Boltzmann ideal gas (the interactions and quantum 
statistics effects are neglected) in the volume V and at temperature T the GCE and CE 
partition functions read: 



oo oo 



z g ^ e xT,v,fi Q ) = y Y 



(X+Z) N + (X-Z) N u Q {N+-N_)/T 



N+=0N-=0 

oo oo 



Y Y Z N+,N-(T,V, fig) = exp (2z cosh [/xq/T]) 

N + =0 V_=0 



(3) 



oo oo 



z c . e .(T,v, Q ) = y y (j ^^^ s (q-i n +- n -]) 

N + =0N_=0 + ' 

oo oo 

= E E Zn + ,nAT,V,Q) 

N+=0N-=0 

= -L f dcj) exp [iQ4> + z (A+ e<* + A_ e*) ] = J (2*), (4) 



where z is a single particle partition function: 



gv_ 

2tt 2 



p 2 <ip e" 



(5) 



g is a degeneracy factor (number of spin states), m - particle mass and \± are auxiliary 
parameters that will be set to unity after calculation of average values. We also labelled the 
number of particles in GCE as (N±). Let us omit the indexes c.e., g.c.e., etc., for partition 
function as the arguments of Z already show to what ensemble it corresponds. The average 
values in both the GCE and CE can be calculated as follows: 



oo oo 
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dz 



(6) 
(7) 



A±=l 



1 e.g. 7r + and 7r mesons 



4 



In thermodynamic limit, V — > oo, and for Q = it gives: 



(N ± ) = z 



1 



{Nl) = z + z\ 

2 \ rs_/ 2 



From the definition of the scaled variance ([I]) it then follows 

_ (Nl)-(N ± f 1 
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FIG. 2: The scaled variances in GCE (horizontal dashed line) and in CE (solid line) for Q = 



Thus for zero system charge in thermodynamic limit the scaled variance in CE is two times 
smaller then in GCE while average particle numbers are the same, see Eqs. (jHJ) and ([9]) left, 
and Eqs. ( flOi) . ( flTi) . One can also see from Fig. [2] that the thermodynamic limit is reached 
very quickly. The scaled variance uf e almost reaches its limiting value at (N±) — z ~ 5-i-10. 

Multiplicity fluctuations for non-zero charge in multi component system with two exactly 
conserved charges, namely electric charge and baryon number, are considered in 9J. The 
relation o; c . e , = oo g _ c _ e _/2 is preserved in multi component system if the number of all positively 
and all negatively charged particles of different species is the same. Large non-zero charge 
Q > leads to additional suppression of u;+ e and the enhancement of ui~ e _, while the relation 
u c.e. < u c.e. < U ^.c.e. holds. Additional baryon charge conservation leads to even stronger 
suppression of the scaled variance in CE comparing to GCE in thermodynamic limit. 
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B. Microcanonical ensemble 



The microcanonical partition function can be easily calculated analytically for the system 
of N noninteracting massless neutral particles if one neglects the effects of quantum statistics. 
This is just iV-times integrated over momentum 5-function 

N ^ nnn / N 



Z N (E,V) 



gv 



iTi- 



p\dpi ... j p 2 N dp N 6\E- ^Pj 

i=i 



1 fgV 



N £3N-1 



(12) 



N\\ix 2 ) (3JV-1)! 

where E - is the energy and V -^volume of the system. One can also generalize Eq. (|12l) to 
the system of charged particles 



J N+,N. 



(E,V,Q) 



1 faV\ N++N - E 3( - N + +N ^- 1 

11 \ = — 8(Q — [N + — N-]) , (13) 



NJ.NJ. 
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and calculate corresponding scaled variances using Eqs. (|T2l) and (|T3l) similarly to (jSJ), ©• 
In thermodynamic limit, V — > oo, and for Q = it gives 



^g.m.c.e. — ^ I 1 



+ . . . 
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+ ... 



(14) 



(iV) ' "V ^~8V^ H52 (Nl) 

where (N) and (A^.) are the average number of particles in GCE. Thus, one can see that 
the scaled variance in thermodynamic limit is 4 and 8 times smaller than in GCE for GMCE 
and MCE correspondingly, see Fig. [31 
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FIG. 3: The scaled variances in the GMCE, left, and in MCE, right [8] 



It means that the thermodynamic equivalence for mean particle number does not apply 



to fluctuations measured in terms of the scaled variance 



(N) ~ (N) c . e . ~ (JV)m.c.e. , V - OO (15) 

Wj.ce. 7^ W c . e . 7^ W m .ce. , V > OO (16) 



see also [lOj for the summary of some limiting values of the scaled variance. Note, that 
average particle numbers in GMCE, MCE and GCE are equivalent in thermodynamic limit 
8] similarly to CE, see (|SJ), (J9J), left. Canonical and microcanonical suppression {§]], 
9] and even microcanonical enhancement [ll[ of average multiplicity (N) is observed for 
very small systems only. Quantitatively, the limiting behavior in the MCE is reached even 
quicker than in CE: for 2 -f- 3 particles if we consider (N) or (N±) and for 3 -j- 4 particles if 
we consider scaled variance see Fig. [3] and [8( . 

The analytic calculations presented above are possible only for Boltzmann statistic in 
CE and for Boltzmann massless particles in MCE. The inclusion of other conserved charges 
and quantum statistic makes the calculations technically very difficult. The simplest way to 



overcome these difficulties is to consider multiplicity distributions in different ensembles 12]. 



III. MULTIPLICITY DISTRIBUTION 



Multiplicity distribution 2 , partition function, different moments, variance and scaled vari- 
ance are closely related, namely: 



P(N) = Z -f , 



{(AN) 2 ) = (N 2 ) - (iV) 5 



(N k ) = ^N k P(N) 



N 



((AAQ 2 ) 
(AO 



Multiplicity distribution P{N) in ideal gas tends to Gaussian Pq(N) for N ^> 1: 



P(N > 1) ~ P G (N) 



1 



exp 



2uj- (N) 



(17) 
(18) 



(19) 



y/2iru}- (N) 

One can easily check this for Eqs. (jlj) and (fT2|) . (TT3T) . see the result in Fig. 0] and 
detailed calculations in CE ]7], MCE and GMCE 8]. 

One can see that multiplicity distributions in different ensembles have the same maximum 
at iV = (N), but different width 3 . As an example in Fig. H]we choose A^ = (N) = 10. One 



2 probability to find some number of particles N if their average number (N) is fixed by external conditions. 

3 in non-relativistic case N = const by definition and P(N) ~ S(N) i.e. it would be a vertical line in Fig.[4j 



7 




can also see that the distributions are smooth and have Gaussian form. Thus, quantitatively, 
iV = 10 is already big enough to consider Gaussian approximation. 

To generalize our formalism for several conserved charges and include quantum statistic, 
let us consider again a gas of Boltzmann particles in the CE for simplicity. The GCE and 
CE partition function are also closely related: 



Z{T,V,iiq) = E eWT E Z N+ , N _(T,V,Q) (20) 

Q=-oo N + ,N- 

The substitution of (j3J in fT2"0"|) transforms it to the identity [301 ]: 

oo 

Z{T,V,fi Q ) = E e Q ^/ T I Q (2z) = exp(2zcosh[/i Q /T]) . 

Q=— oo 

After the replacement e^ MQ//T = e( N+ ~ N ~^ Q / T and z± = 2e ±AtQ//T one obtains: 

oo 7V+ iV_ 

z(t,t^ q ) es z = E E ^^iW-[JV + -JV-]) 

Q=-oo7V+,7V_ + ' 
Q=-oo"'~' 7r 

= E / i ?- lQ * m = E z « • < 21 > 

where 2(0) is the GCE partition function with replaced chemical potential Hq/T — > hq/T + 
i<p. Similarly to ffTTj) . the probability of finding the GCE system with the particular net- 



charge Q equals the following [121 ]: 



P(Q) 



z 



2tt 



(22) 



The probability to find the number of positively charged particles N + that is exactly equal 
to iV in the GCE system with net-charge equal to Q is as follows [12I ]: 

1 N + N- 

p (N,Q) = ^ E %h%h6{Q-[N + -N-])6(N-N + ]) 



Z ^ N + \ NJ. 

N+,N- + 
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N\(N-Q)\ 



, (23) 



where Z( 



exp \z (^ e t J -Q/T+i<f>+i<i>N _|_ e mq/t Finally, the particle number distribu- 

12]: 



tion in CE can be found as a ratio of the distributions and 

p(iv,g) 



calculated in GCE 



P(tf|Q) 



One can easily check that 



(JV{> 



-)-/ c.e. 



^iV fc P(iV|Q) 

iV 



V 7V „iV- 



iV! (N-Q)\ 



(24) 



(25) 



The Eq. ( I24p is very important, because it allows to calculate a value in CE using the values 
calculated in GCE. It also allows for generalization to quantum statistic and taking into 
account several exactly conserved charges, energy conservation, resonance decay, etc. To 
do this one just need to take corresponding GCE partition function and multiply it by the 
Fourier representations of the relevant delta functions 121 ]: 



P(Q j ) = % 



J oi 1 



(26) 



where j runs over all conserved quantities. Repeated upper and lower indexes j imply 
summation over j. The function Z(4>) also changes if we include different particle species 
and quantum statistic: 



Z(4>j) = exp 



(27) 



where the single particle partition function of particle specie I is given by: 



gjv 

(2vr) 3 



d 3 p In 



1 ± e~( £i ~ w )/ T e^^j 



±1 



EE V^^). 



(28) 



We introduced here particle Vs charges q\ = qi = (qi, bi, sj, ...) that corresponds to the charges 
conserved in the system. We also introduced the degeneracy factor gi = (2J; + 1), internal 
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angular momentum J;, mass mi, and energy e\ = \Jp 2 + mf, the chemical potential vector 
/ iJ = (/Wq, A*5---)> an d particle Z's chemical potential \ii = qjfij. V is the system volume, 
and T it's temperature. The summation Y2i includes also anti-particles, for which qj — > —q^. 
The upper sign in the Eq. (I28p denotes Fermi-Dirac statistics, while the lower is used for 
Bose-Einstein statistics. The Boltzmann approximation is obtained from (1251) as a first term 
of the series expansion for e~ ( - £l ~'^ l ^ T <C 1. 

The real calculations of (126]) can be performed only in the limit V — > oo. Then the main 
contribution to the integral in (1261) comes from a small region around the origin. Thus it is 
possible to make the Taylor expansion of Yli ipi an d leave only the first two terms. Similar 
saddle point expansion was intensively used for partition function itself [l?]], {21], 22], 23] 
while the relations between partition function, multiplicity distribution, and scaled variance 
was obtained only in [l^ . 

It was shown for GMCE i n | l3l| and for the most general case of MCE with arbitrary num- 
ber of conserved charges in [12| that the variance is proportional to the ratio of correlation 
matrix determinants: 



((AiV) 2 ) 



V 



det 


A 


det 





(29) 



where the elements of the correlation matrixes can be found as follows: 



.4 



d 2 \ogZ{ 



Ai 



d 2 \ogZ( 



0=0 



(30) 



4>,(f> N =o 



The difference between A and A is that in the latter case i and j run over iV also 
the scaled variance is a ratio of (|30|) to the mean multiplicity: 

. d log Z (0 j,(f) N 



12| . Then 



(N) 



(31) 



<f>,<t>N=0 

The above method is very powerful. Nevertheless it fails in the case of Bose condensation, 
because scaled variance in GCE then goes to infinity [ijj] and multiplicity distribution has 
infinite width. All matrix elements (13"01 and higher derivatives of Z(<f)) tends to infinity 
in GCE jl2]. However, exact charge and energy conservation suppress even these infinite 
fluctuations [l4], [l5]. The very special selection of events is need to see these infinite 
fluctuations in MCE. This is proposed as the signal of possible 7r-meson condensation in 
p + p collisions Jj]] , • 

The further improvement is possible if one consider average multiplicity and fluctuations 
at different momentum levels. This approach is called the microcorrelator method [sj]. It 
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analogous to the above approach [12|], but additionally allows to consider correlations be- 
tween different momentum levels. The full hadron gas in the next section is considered using 
microcorrelator method [l7], j^ol ]. 

IV. HADRON GAS 

Let us consider the fluctuations in the ideal relativistic gas with different types of hadrons 
in the MCE with exactly fixed the global electric (Q), baryon (B), and strange (S) charges 
of the statistical system. The system of non-interacting Bose or Fermi particles of species 
i can be characterized by the occupation numbers n p ^ of single quantum states labelled by 
momenta p. The occupation numbers run over n Pi i = 0,1 for fermions and n Pt i = w x 1,2,... 



for bosons. The GCE average values and fluctuations of n p $ equal the following 



29|: 



( n v,i) = — r/ ; 9 X \ — I ' ( 32 ) 

a/p 2 + mf - Hi) /T 



cxp 



7i 



2 ~ ( An l,i) = (i n P,i ~ ( n P,i)f) = ( n P,i) ( x + 7i<np,i» • ( 33 ) 



v p,* 



In Eq. (|32|) . T is the system temperature, m ; is the mass of i-th particle species, 7« corre- 
sponds to different statistics (+1 and —1 for Bose and Fermi, respectively, and 7« = gives 
the Boltzmann approximation), and chemical potential /ij equals: 

A*i = Qi + h + Si /j s , (34) 

where qi, bi, Si are the electric charge, baryon number and strangeness of particle of specie 
i, respectively, while /xq, fi B , fx s are the corresponding chemical potentials which regulate 
the average values of these global conserved charges in the GCE. 

The average number of particles of species i, the number of positively and negatively 
charged particles are equal: 

(m = = £2 / p 2d p m . ( N +) = E w - ( N -) = E w ' 

p i,qi>0 i,qi<0 

(35) 

where gi is the degeneracy factor of particle of species i. A sum of the momentum states 
means the momentum integral, which holds in the thermodynamic limit V — > 00. 

Particle number fluctuations and correlations can be calculated in all ensembles using the 
microscopic correlator method. 

(ANi ANj )... =J2 ( A n p ,i An k>j )... , (36) 

p.k 
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where ( )... means GCE, CE, or MCE microscopic correlator. The scaled variances of 
negatively and positively charged particles read: 



LO 



((AiV-) 2 ) 
<AL> 



LO 



((AiV + ) 2 ) 



(37) 



where 



((AN.) 2 ) = (ANiANj) , ((AN + ) 2 ) = £ (A^AN,) . (38) 

qi<Q,qj<0 qi>Q,qj>0 



The microscopic correlator in the GCE reads: 



(An p>i An kJ ) 



(39) 



where v 2 i is given by Eq. fl33|) . This gives a possibility to calculate the fluctuations of 
different observables in the GCE. Note that only particles of the same species, i = j, and 
on the same level, p = k, do correlate in the GCE. Thus, Eq. ( 1391) is equivalent to Eq. ( I3B1 : 
only the Bose and Fermi effects for the fluctuations of identical particles on the same level 
are relevant in the GCE. 

The MCE microscopic correlator is as follows [l^], 2( 



V p,i V k,j 



(An Pi jAn feJ ) m . c . e . = v p x 5 pk TW 1 [ 1i1j M Qq + b i b j M bb + SiSjM s , 

+ (liSj + QjSi) M qs - (q,ib,j + qjbi) M qb - (biSj + bjSi) M bs 

+ e pi e kj M ee - (qitpj + qjt ki ) M qe + (bie pj + bjt ki ) M be - (si€ pj + Sje ki ) M s , 

where \A\ is the determinant and My are the minors of the following matrix: 

(A(q 2 ) A(bq) A(sq) A(eq)\ 



(40) 



A 



A(qb) A(b 2 ) A(sb) A(eb) 
A(qs) A(bs) A(s 2 ) A(es) 
\A(qe) A(be) A(se) A(e 2 ) J 



(41) 



v 



p,k ' 



with the elements, A(g ) = ^, fc q k v P ,k > A (Q b ) = E P ,fc^ fo ^p,fc . A (<? e ) = E P , fc Qke pk 
etc. The sum, k , means integration over momentum p, and summation over all hadron- 
resonance species k contained in the model. Note that the presence of MCE terms containing 
single particle energies, e pi - 



lp 2 + m 2 , in the last line of Eq. po!) is a consequence of exact 
energy conservation. In the CE, only charges are conserved exactly, thus the terms of the 
last line in Eq. fl4"0l are absent, and A in Eq. (14"T|) becomes the 3x3 matrix (see Ref. jl?!]). 
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V. EFFECT OF RESONANCE DECAYS 

The average number of 2-particles in the presence of primary particles N* and different 
resonance types R is the following: 

(Ni) = (N*) + J2( N R)J2 b r<r = W> ( 42 ) 

R r R 

The summation YIr, runs over all types of resonances. The (. . .) and (. . .) R correspond 
to the GCE averaging, and that over resonance decay channels. Resonance decay has a 
probabilistic character. This itself causes the particle number fluctuations in the final state. 



In the GCE the final state correlators can be calculated as 24j : 



(ANiANj) = (AN* AN*} + ^ [(AN%) (m) R (n 3 ) R + (N R ) (A^An,-)*] , (43) 

R 

where is the branching ratio of the r-th branch, nf r is the number of 2-th particles 
produced in that decay mode, and r runs over all branches with the requirement Y2 r ^ — 1 
and (Arij Arij) R = ^ r bfnf r nf r — (ni) R (rij) R . Note that different branches are defined in 
a way that final states with only stable (with respect to strong and electromagnetic decays) 
hadrons are counted. 

All primary particles and resonances become to correlate in the presence of exact charge 



conservation laws. Thus for the MCE correlators we obtain a new result 17| : 
(ANiANj)^. = (AN*AN*) m . c . e . + ^(N R ) (An t An,} R + ^(AiV* AN R ) m . c . e . ( nj ) R 

R R 

+ ^(AiV* AN R ) m . c . e . (m) R + J2^ ANr AN R')m.c. e . (ni) R (n,)^ . (44) 

R R,R> 

Additional terms in Eq. ( j44l compared to Eq. (f4*3l) are due to the correlations induced by 
exact charge conservations in the MCE. The Eq. (jHJ) remains valid in the CE too with 
(• • )m.c.e. replaced by (. . .) c . e ., the difference between them appears only when one specifies 
the microscopic correlators (f40|) of the MCE or CE. 

VI. SCALED VARIANCES ALONG THE CHEMICAL FREEZE-OUT LINE 

Mean hadron multiplicities in heavy ion collisions at high energies can be approximately 
fitted by the GCE hadron-resonance gas model. The fit parameters are temperature T, 
chemical potentials (/is, fj,s, £*q), and strangeness suppression factor 75-, which allows for 
non-equilibrium strange hadron yields. There are several programs designed for the analysis 
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of particle multiplicities in relativistic heaw-ion collisions within the hadron-resonance gas 
model, see e.g., SHARE [25(] , THERMUS 26j and THERMINATOR [27|. In this paper an 



26j is used. 



extended version of the THERMUS thermal model framework 

For the chemical freeze-out condition we choose the average energy per particle 

n n 

(E) j (N) = IGeV [28]. Using the standard parametrization [5] we obtain the T — /ie 
freeze-out line for central A+A collisions (see Fig. [5]). The center of mass nucleon-nucleon 
energies, \/Snn, marked in the figures below correspond to the beam energies at SIS (2A 
GeV), AGS (11. 6A GeV), SPS (20A, 30A, 40A, 80A, and 158A GeV), colliding energies at 
RHIC (v 7 ^ = 62.4 GeV, 130 GeV and 200 GeV) and LHC (v^vv = 5500 GeV). 
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FIG. 6: The scaled variances for negatively and positively charged particles, both primordial and 
final, along the chemical freeze-out line for central Pb+Pb (Au+Au) collisions. Different lines 
present the GCE, CE, and MCE results. Symbols at the CE and MCE lines for the final particles 
correspond to the specific collision energies. The arrows show the effect of resonance decays 2o| . 
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Figure [6] show the prediction for the scaled variances for negatively and positively charged 
particles as a function of ^snn- 

The prediction can be compared with the preliminary NA49 data on Pb+Pb collisions 
at 20A-158A GeV [3] using the following approximate formula: 



(45) 



where u i7T refers to an ideal detector with full 47r-acceptance and a>* c is the scaled vari- 
ance measured by a real detector with a limited acceptance), q is the ratio between mean 
multiplicities of accepted particles and all hadrons. In the limit of a very 'bad' (or 'small') 
detector, q — > 0, all scaled variances approach linearly to 1, i.e., this would lead to the 
Piossonian distributions for detected particles. However, we find a strong qualitative differ- 
ence between the predictions of the statistical model valid for any freeze-out conditions and 
experimental acceptances: the CE and MCE correspond to c e < u;^ e < 1, and the GCE 
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FIG. 7: The scaled variances for negative (left) and positive (right) hadrons along the chemical 
freeze-out line for central Pb+Pb collisions at the SPS energies. The corresponding T and \xb 
values at different SPS collision energies are presented in Fig. [SJ Different lines show the GCE, 
CE, and MCE results calculated with the NA49 experimental acceptance 



From Fig. [7] it follows that the NA49 data for uj^ extracted from the most central Pb+Pb 
collisions at all SPS energies are close to the results of the hadron-resonance gas statistical 
model within the MCE. The data reveal even stronger suppression of the particle number 
fluctuations. A possible reason of this is an uncertainty in the determination of the detector 
acceptance and an additional suppression due to momentum conservation and the excluded 
volume effects in the hadron-resonance gas. 
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In order to allow for a detailed comparison of the distributions the ratio of the data and 
the model distributions to the Poisson one is presented in Fig. [HJ The convex shape of 
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FIG. 8: The ratio of the multiplicity distributions to Poisson ones for negatively charged hadrons 
produced in central (1%) Pb+Pb collisions at 20A GeV, 30A GeV, 40A GeV, 80A GeV, and 158A 
GeV (from left to right) in the NA49 acceptance The preliminary experimental data (solid 
points) of NA49 [19] are compared with the prediction of the hadron-resonance gas model obtained 
within different statistical ensembles, the GCE (dotted lines), the CE (dashed-dotted lines), and 



the MCE (solid lines) 



20] 



the data reflects the fact that the measured distribution is significantly narrower than the 
Poisson one. This suppression of fluctuations is observed at all five SPS energies and it is 
consistent with the results for the scaled variance shown and discussed previously. The GCE 
hadron-resonance gas results are broader than the corresponding Poisson distribution. The 
ratio has a concave shape. An introduction of the quantum number conservation laws (the 
CE results) leads to the convex shape and significantly improves agreement with the data. 
Further improvement of the agreement is obtained by the additional introduction of the 
energy conservation law (the MCE results). The measured spectra surprisingly well agree 
with the MCE predictions [20(. 



VII. SUMMARY 



We have found that scaled variances are different in different statistical ensembles. For 
relativistic one component Boltzmann gas with zero charge in thermodynamic limit we 
analytically obtained rather interesting limiting values: w 9 . c . e . = 1, u c .e. = 1/2, <*>j gm . c . e (rn = 
0) = 1/4 and oJ m , c , e ,{m = 0) = 1/8. We also found an analytical method to account for 
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resonance decays. The formalism that allows to consider any number of conserved charges 
and also energy conservation in full hadron-resonance gas was developed. 

The experimental data allows to exclude GCE for scaled variance. They show reason- 
able agreement with CE and surprisingly well agree with the expectations for the MCE. 
Thus the predicted suppression of the multiplicity fluctuations in relativistic gases in the 
thermodynamic limit due to conservation laws do exist. 
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